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Abstract 



(N 

^ ■ We study the pricing problem for a European call option when the volatility of the underlying 

oo . asset is random and follows the exponential Ornstein-Uhlenbeck model. The random diffusion 

in 

^N \ model proposed is a two-dimensional market process that takes a log-Brownian motion to describe 

CD ' price dynamics and an Ornstein-Uhlenbeck subordinated process describing the randomness of the 

OO ■ 

^D ■ log- volatility. We derive an approximate option price that is valid when (i) the fluctuations of the 



volatility are larger than its normal level, (ii) the volatility presents a slow driving force toward 



H ' its normal level and, finally, (iii) the market price of risk is a linear function of the log- volatility. 

We study the resulting European call price and its implied volatility for a range of parameters 
consistent with daily Dow Jones Index data. 
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I. INTRODUCTION 



The picture that asset prices follow a simple diffusion process was first proposed by 
Bachelier in 1900 and his main motivation was precisely to provide a fair price for the 
European call option. The buyer of the call option has the right, but not the obligation, to 
buy the underlying asset from the seller of the option at a certain time (the expiration date) 
for a certain prestablished price (the strike price) pj. 

After the introduction in 1959 of the geometric Brownian motion as a more refined market 
model for which prices cannot be negative, the volatility can be viewed as the diffusion 
coefficient of this random walk [2] . This simple assumption - the volatility being a constant 
- lies at the heart of Black-Scholes (BS) option pricing method 3|. Within the BS theory, 
the most direct technique constructs an equivalent martingale measure for the underlying 
asset process. This change of the probability measure guarantees a fair price for the option. 
The price of the European call is readily obtained: it simply consists of an average over the 
final payoff evaluated under the martingale measure, and appropriate discounting. 

However, and especially after the 1987 crash, the geometric Brownian motion model and 
the BS formula were unable to reproduce the option price data of real markets. Several 
studies have collected empirical option prices in order to derive their implied volatilities (the 
volatility that agrees with the BS formula). These tests conclude that the implied volatility 
is not constant and seems to be well-fitted to a U-shaped function of the ratio between the 
price of the underlying asset and the strike price of the option (the so-called "moneyness" ) . 
This phenomenon, known as the smile or smirk effect, the latter term due to the fact that 
the U-shape is usually asymmetric, shows the inadequacy of the geometric Brownian model. 
Nevertheless, the continuous-time framework provides several alternative models specially 
designed to explain, at least qualitatively, this effect. Among them, we highlight on the 
Stochastic Volatility (SV) models. These are two-dimensional diffusion processes in which 
one dimension describes the asset price dynamics and the second one governs the volatility 
evolution (see e.g. Refs. [J, bl. Id, U, M, Id, M llll, M)- 

The SV class of models have also been studied to check whether they are also able to 
describe the dynamics of the asset itself. In this direction there are several interest ing 



works in the literature focusing on t 
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lis and others similar issues (see e.g. Refs. 13 
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251]). It is indeed worth noticing that the theoretical 



framework of SV modeling has many common points with the research on random diffusion 
whose aim is to describe the dynamics of particles in random media which can explain a 
large variety of phenomena in statistical physics and condensed matter 26|]. Going back 



to finance, there is a wide consensus t 
financial markets 
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lat volatility plays a key role in the dynamics of 



33| . Among the most relevant statistical properties 



of the financial markets, volatility seems to be responsible for the observed clustering in 
price changes. That is, large fluctuations are commonly followed by other large fluctuations 



and similarly for small changes 



27(1 . Another feature is that, in clear contrast with price 



changes which show negligible autocorrelations, volatility autocorrelation is still significant 



for time lags longer than one year ll|, ll3|, |27|, 



28 



34| . Additionally, there exists the 



so-called leverage effect, i.e., much shorter (few weeks) negative cross- correlation between 



current price change and future volatility [Ij, |27|, |31|, l32j. Finally, several authors have 
argued that a g ood approximation for the volatility distribution can be given by the log- 



normal 



20 



27 



30|. Another possible candidate to such distribution is provided by the 



so-called inverse Gaussian distribution [33]. 

All these observations have led us to consider the exponential Ornstein-Uhlenbeck (ex- 
pOU) model, since we have seen that, among the simplest and classic SV models, the expOU 
model is able : (i) to describe simultaneously the observed long-range memory in volatility 



and the short one in leverage 
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16 



22| , (ii) to provide a consistent stationary distribution 



for the volatility with data 12, l20|, l23|] , (iii) it shows the same mean first-passage time profiles 
for the volatility as those of empirical daily data [2^ , and finally (iv) it fairly reproduces the 
realized volatility having some degree of predictability in future return changes 23|]. Our 
aim in this research is to take advantage of all this knowledge to study the European option 
price. We shall propose an approximate price for the option valid for realistic parameters 
that guarantees the empirical observations just mentioned. We therefore study the influence 
of these empirical observations in the price by designing an approximation procedure that 
provides an alternative price which is appropriate for a different range of parameters to those 
already published in Refs. (Sj (see also the discussion at the end of Sect. III). Moreover, our 
present study goes much further than the simulation analysis already performed by one of 

us y. 



Before closing this introductory section, we should remark that the approximation proce- 
dure we use has been specially designed to tackle this kind of financial problems. However, 



the method is somewhat related to the Born-Oppenheimer approximation which, after its in- 
troduction in 1927 for the quantum mechanical treatment of molecules, has been widely used 
for a great variety of physical applications. In statistical mechanics the Born-Oppenheimer 
approximation is specially suited for the so-called "adiabatic elimination of fast variables" 



(in our case the volatility) 35j. A simplified version of that approximation has been recently 
used to address some financial problems [l5|, [l9| . 

The paper is divided into five sections. Section [TTl presents the expOU model and shows 
its main statistical properties. In Section UTTl we construct the equivalent martingale density 
and afterward derive the European call option price. Section IIVI studies the price derived 
and compare it with the classic BS price. Conclusions are drawn in Section |V] and some 
technical details are left to the Appendices. 

II. THE MARKET MODEL: A SUMMARY 

As most SV models the correlated exponential Ornstein-Uhlenbeck (expOU) stochastic 
volatility model is a special kind of a two-dimensional diffusion process. We have thoroughly 
studied the statistical properties of the model in previous works with a quite satisfactory 



agreement with empirical data of real markets 12, l2l|, l22|, |23|, l2J] . The model proposed by 



Scott [7| in 1987 is furnished by a pair of Ito stochastic differential equations: 

^ = /irft + me^Wrfl^i(t) (1) 

S{t) 

dY{t) = -aY{t)dt + kdW2{t), (2) 

where S{t) is a financial price or the value of an index. The parameters a, m, and k are 
positive and nonrandom quantities and dWi(t) = C,i{t)dt {i = 1,2) are correlated Wiener 
processes, i.e., C,i{t) are zero- mean Gaussian white noise processes with cross correlations 
given by 

m)ut')) = Pv5it-t'), (3) 

where Ai = 1, Pij = p {i j^ j, —I < p < 1). We remark that while the original Scott 
model [7| had independent Brownian motions, it is important to allow for non-zero pij to 
explain some empirical observations [sl, |9| . 
From Eq. ([2]) we see that 



Y{t) = Foe""^*-*"^ + k f' e-"^'-'UW2{s), (4) 



where we assume that the volatihty process Y{t) starts at certain initial time to (which 
can be set equal to 0) with a known value Y{to) = Yq. The process Y{t) is Gaussian with 
conditional first moment and variance given by 



2a 



E[Y{t)\Yo] = roe-"(*-*°) and Var[F(t)|yo] = — (l - e" 
In the stationary limit, {t — to) -^ oOj we have 



2a{t-to) 



(5) 



E[Y{t)] = and Var[r(t)] 
In terms of the process Y{t), the volatility is given by 



2a' 



(6) 



a(t) = me^W. 
Hence, the conditional probability density function (pdf ) for the volatility is 



(7) 



p(o-,t|ao,0) 



and the stationary probability density reads 

. . 1 



exp 



[\n{a/m) — e "*ln(o"o/m)]" 
2/32(1 - e-2"*) 



(8) 



Pst[(^) 



where 



ay/2^ 



exp 



{-ln\a/m)/2f3^} 



2a' 



(9) 



(10) 



This lognormal curve for the volatility is consistent with real data as it has been reported 



in several studies for different time lags 12, 



20 



23 



21, 



3Q| 



Other interesting statistical properties of the model refer to the returns process given by 

dR{t) = dS{t)/S{t) - E[dS{t)/S{t)]. 



We thus have the squared returns autocorrelation [l2| 

exp[4/?2e-°^] 



Corr 



dR{ty,dR{t + T) 



3e4/3' - 1 



(11) 



and the leverage effect (or return-volatility asymmetric correlation) 12| 
E [dR{t)dR{t + r)2] 2pk 



C{r) 



m 



■ exp [-ar + 2/5^(6-"" - 3/4)} iJ(r), 



(12) 



E [dRitfY 
where H{t) is the Heaviside step function (if (r) = l when r > 0, and H{t)=0 when r < 0) 



Both the volatihty and the squared returns process (fTTj) show a decay in their autocorre- 
lation described by the desired cascade of exponentials which goes from very fast time scales 
to the slowest ones 



3J|. From Eq. (ITTj) . one can obtain 



Corr 
which goes as 

Corr 



dR{ty,dR{t + rf = .,^, . ^ 



1 ^ (4/? 



2\n 



{3e^0^-l)^^ n\ 



— nOLT 



dR{t)\ dR{t + r)2] ^ ;_ e-"^ for ar » 1. (13) 



For a small enough, of the order of 10^'^days^ , this may appear as a long range memory 
up to few trading years [12] ■ 

The model is similarly able to explain the rather swift leverage effect -up to few weeks- 



which is well observed in actual financial markets. From Eq. (IT2|) . one can now obtain 12| 



C(t) ~ ^ exp(f3y2)e-'''^H(T) for ar < 1, (14) 

m 



and 



£(r) ~ ^exp(-3/?V2)e-""i7(r) for ar > 1. (15) 

m 

Note that the leverage amplitude in the latter regime is smaller than that over the shorter 

range (cf. Eqs. ( TT^ and (IT5l)). Hence, the leverage asymptotics given by Eq. ( IT5|) is negligible 

compared to the one provided by Eq. 01 4p . This result is consistent with real data showing 

that the leverage only persists up until few weeks and thus has a smaller range than the 

volatility autocorrelation. In this case, the model provides an exponential decay consistent 

with data if one takes characteristic time scale (given by the k"^ if we assume k'^ > a) 

between 10 — 100 days and a moderate negative correlation between the two Wiener input 

noises p ~ —0.5 [12| . 

We remark that the range of parameters mentioned above is outside the values taken in 

Ref. [8|. In that case, Fouque et al assume that /5^ ~ 1 (so that fc^ ~ a), thus implying that 

leverage and volatility autocorrelation have a very similar decay with a characteristic time 

of the same order (cf. Eqs. flTT|) - flT2|) and flT3l) - flT4l) ). Herein we consider a broader range 



for /3^ 12| in order to get a leverage decay faster than the volatility autocorrelation k'^ ^ a. 
Empirical observations have led us to consider the range where /?^ is much larger than one. 
In addition, Fouque et al also take a such as to provide a fast reversion to the normal level 







. Even 



of volatility with a characteristic time scale 1/a between one day and few weeks 

from this perspective their framework is far away from the one we here propose since we are 

assuming that reverting force is quite slow with l/a of the order of one or more years. 

III. THE OPTION PRICING PROBLEM 

In a previous paper pJi] we have addressed the question of the distribution of returns 
characterized by either the return probability density function (pdf) or by its characteristic 
function with a rather satisfactory fit for 20 day price return changes. We have there derived 
an approximate solution for the return pdf in terms of an expansion of Hermite polynomials 
with prefactors related to the skewness and the kurtosis of the expOU model. We want now 
to apply a similar analysis to the option pricing problem illustrated by the European call 
C{S,T). The payoff of this contract, that is, its final condition, is 

C(5,0) = max[5(T)-ir,0]. (16) 

where S{T) is the underlying asset price at expiration date T and K is the strike price. 

A. The equivalent martingale measure 

Following the standard approach to option pricing [l|, l9| , we pass to an equivalent mar- 
tingale (or risk-neutral) pricing measure P*, under which the discounted price of traded 
securities are martingales. In particular, assuming a constant risk-free interest rate r: 



rt 



E*[S{t)\S] = Se 

where 5* = 5'(0). Under P*, the dynamics ([I])-© of the expOU model are given by the 
following set of stochastic differential equations 

^ = rdt + me^^''>dW*{t) (17) 

S{t) 

dY{t) = -aY{t)dt-kA{X,Y,t)dt + kdW;{t), (18) 

where {W^, W2) are P*-Brownian motions with the same correlation structure as the original 
Brownian motions {Wi, W2). 

In writing Eqs. flT71) - (TT8ll one has to include an arbitrary function A(-) called the log- 



volatility's market price of risk 



flig 



9|. This function quantifies the risk aversion sensitivity 



of the trader to the volatihty uncertainty. Note that if there were no volatihty fluctuations 
[i.e., k = 0) the function A would vanish. 

It does not seem possible to obtain the exact solution to this option problem and we will 
therefore search for approximate expressions but using different methods than those already 
proposed in the literature [Sl, |9(]. In order to specify the risk aversion function A(X,Y,t) 
we first suppose that it depends solely on the volatility level and not on the return or the 
time. That is, A{X,Y,t) = A{Y). Moreover, we further assume that this dependence is 
weak which amounts to A{Y) being approximated by a linear function of Y: 

A{Y)^Ao + A^Y. (19) 

We now define a new volatility-driving process 

Z = Y + ^, (20) 

a 

and the "log return" X{t) by 

X{t) = \n[S{t)/S]. (21) 

Then Eqs. fITTj) and flTSl) turn into 



dX{t) 



- -m2e2^(*) 



dt + me^WdW^*(t) (22) 



2 
dZ{t) = -aZ{t)dt + kdW;{t), (23) 

where 

fh = mexp{—k Aq/ a) , a = a + kAi. (24) 

We analyze the consequences of the introduction of a risk aversion function in the equiva- 
lent martingale process. Note first that the parameter a which quantifies the memory range 
of the volatility (cf. Eq. ( IT3l) with a replaced by a) increases when Ai > 0. Secondly, when 
Aq > the normal level of volatility shifts to lower values and the stationary volatility pdf 
has a smaller standard deviation (cf. Eqs. Qj, flTOj) and flM|) ). 

B. Approximate solution of the Fokker-Planck equation 

Let us denote by P2ix, z, t|0, ^o) the joint density of (Xj, Zt) under the equivalent martin- 
gale measure. This density obeys the following Fokker-Planck equation 



with initial condition 

P2{x, z, 0\0, zo) =6{x)6{z- zo). (26) 

Even when the risk aversion is a hnear function of the log-volatihty, as is our case, we 
may only derive an approximate solution valid under particular values of the parameters of 



the model for a given market [l2l]. Herein we will treat the case in which the "vol of vol" k 

is much greater than the normal level of volatility rh. In other words, we will assume that 

the parameter 

A = - » 1 (27) 

m 

is large. Certainly this is at least the case of the Dow Jones Index daily data for which 

A ~ 102 [12]. 

We define a dimensionless time t' and two scaling variables u and v by 

t' = kH, u = Xx, v = Xz. (28) 

The FPE ([25]) then reads 

dt' dv^^^'^2dv^ X\m^ 2 ) du 



and 



P2{u,v,Q\{),vq) =5{u)5{v -vq), (30) 

where 

Vq = Xzo, and u = 1/(2/5^) = a/k'^. (31) 

We can write a more convenient equation in terms of the characteristic function defined 

as 

/OO POO 

e^'^'^'du e"^^''p2{u,v,t'\{],VQ)dv. (32) 

-OO J — OO 

This double Fourier transform turns the initial- value problem (1291) - (1301) into 
d d 1 

— V5(u;i, UJ2, t') = -UUJ2-K—'^{^1, ^2, t') - -X^UjIlp{uJi, U)2, t') 

Z IT 

-—UiipiuJi, UJ2 - 2ilX, t') + -—^UJiip{uJi,UJ2, t') 

--oolifi {uji, UJ2 - 2i/X, t') - XpuJiUJ2(p {uJi, UJ2 - i/X, t') , (33) 



with 



<^(a;i,a;2,0|0,t;o) = e^"^'^°. (34) 

Since the payoff for the European option depends on the price but not on the volatihty, 
we only need to know the marginal characteristic function of the (martingale) return X{t). 
This would imply solving Eq. (155]1 when 072 = and assume Ui = uj/\ (see Eq. (1251) ). 
Unfortunately, we cannot proceed in such a direct way and have to solve the two-dimensional 
problem when UJ2 is small and A is large. This is done in Appendix \M where we prove the 
following approximate expression for the marginal characteristic function of the return: 



¥?(cj/A, UJ2 = 0, t'\0, vo) = exp [-C{iu/\, t') + O (l/A^) } . 



(35) 



where C{uj/X, t') is given in Eq. flA12p which, after recovering the original variables t and zq 
(cf. Eq. f l28|) ) and some reshuffling of terms, reads 



m^t 



C(a;/A, at) = i^i{t)uj + '^^^uj"^ + §{t, zo)uj'^ - ipq{t, zq)uo^ - K{t)uj^ + 0(1/A^). (36) 



where 



/i(t) = rt m^t. 



d{t,z^) 
^(t, ^0) 

K{t) = 



^0 ,^_^-., 



AV 
1 

A3^ 



2AV3 



at - (1 - 
at + -(l-i 



-at 



-2at 



A3zy2 
2(1- 



ate""* - (1 - e""* 



(37) 
(38) 
(39) 



-at 



+ 



4, ,3 



2A4z/ 



at 



2(1 



-at 



ate 



-at 



(40) 



The approximate expression for the marginal characteristic function of the return is thus 
obtained by taking first terms of Taylor expansion of Eq. ( !35l) . We get 

(p{u/X, at) = exp< — iujfi{t) + m^ujH/2 > 1 — ^{t, zo)uj'^ + ip<;{t, Zq)uj^ 

+ (K{t) + ^(t, zof/2) uj^ + 0(1/A^)] . (41) 

The inverse Fourier transform of this expression yields the following approximate solution 
for the pdf of the (martingale) return: 



p{x,t\zo) 



\/2iTmH 



exp 



(x — /i)' 
2mH 



d f X — fi 

-"2 



2mH 



V2mH, 



/ x-/i \ K + ^'^12 / x-fi Y 



(42) 



10 



where fJ.,^,<^, and k are given by Eqs. (137|1 - (H0|) . and Hn{-) are the Hermite polynomials: 



X 



"g-(«^)^-™a:^^2; 



TT 



.e-'y^-'Hn{w/2a). 



(43) 



-2za)"a 

As shown in Eqs. (!24|) . we note that both the rate of mean-reversion of the log- volatility 
given by a and the normal level of volatility m are rescaled depending on the own risk 
aversion of the agent. We also recall that the dependence on the original log- volatility 
yo = ln(cro/m) is provided by (cf. Eq. (l20l) ) 



kAo 
zo = yQ + -^. 
a 



(44) 



Moreover, a further simplification can be achieved by averaging over the initial volatility 
with the zero- mean stationary (and Gaussian) distribution of the process Z (cf. Eqs. ( I23l) ). 
We can use same structure as given by Eq. (142!) but replace some of the parameters involved 
by 






A3z/2 



at 



1 



K{t) + ^l9{t, Zof ^ k{t) 



1-e 
1 



-at 



2AV3 



at - (l - 



e ^']+p' 



at 



2(1 



e-^ +ate~^' 



(45) 
(46) 

i[ 

(47) 



In Fig. [T] we represent the approximate expression of the equivalent martingale measure 
as given by Eq. fj42l) . The initial log- volatility is zq = and the linear market price of 
risk is characterized by Aq = Ai = 0.001 (cf. Eq. (TT9|) ). We plot a set of distributions 
where we provide three different values of each parameter (fc, m, p, and a) defining the 
expOU model. We observe that the vol-of-vol k mainly modifies the positive wing of the 
distribution. The normal level of volatility m broadens the probability distribution with 
no difference between the two tails. The correlation between Wiener noises p provide the 
observed negative skewness only if p is negative so that this term should be taken into account 
if we want to include this effect to the corresponding option price. And finally, the long- 
range memory parameter a of the reverting force has little effect on the distribution profile. 
Hence overestimating (or underestimating) this quantity does not have great consequences 
in providing a good approximation of the risk-neutral distribution. 
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FIG. 1: Risk- neutral return density (j42l) for t = 20 days and with terms provided by Eqs. ([37])- 
lOl) when zq = and assuming Aq = 10"^ and Ai = 10^^ (cf. Eq. (I19p ). We depart from the 



parameters m = 10 ^ day 



-1/2 



a 



8 X 10"^day"\ 



P 



-0.4 and k = 0.11 day ^/^ and shghtly 



modify them in each of these plots . 



C. The European Call price 



Our main goal is to study the European call option price. Once we have obtained our 
approximate solution for the equivalent martingale measure, the price can be computed in 
terms of the expected payoff (I16p under our equivalent martingale measure [l| . That is: 



C{S,T,zo) 



-rT 



E 



-rT 



max(5e^(^) - K, 0) X(0) = 0, Z(0) 

/oo 
max(S'e^ — K, 0)p(x, T\zQ)dx, 
-oo 



Zq 



(48) 



where the price at the expiration date T is provided by the return path according to the 
relation S{T) = S exp{X{T)) and K is the strike price. As can be observed from the 
expansion (j^2l) . the computation of the approximate option price implies evaluating four 
integrals. The first one, Cbs, corresponding to the classic Black-Scholes price {i.e., when 



12 



the underlying process has a constant volatihty given by m): 

Cbs{S,T) = SN{di) - Ke~'-^N{d2), (49) 

where 

\nS/K + (r + my2)T , , \nS/K + (r - my2)T , , 

di = , and d^ = ._ (50) 

and N[d) is the normal distribution 

N{d) = -^ I'' e-'^'/^dx. (51) 

\/2,TT J —CO 



. 27r "'-oo 

In the Appendix [B] we evaluate the rest of terms. Summing them up, our approximate 
price for the European call option when underlying follows an expOU stochastic volatility 
process reads 



C{S, T, zo) =Cbs{S,T)+{^ + p<; + >,+ -) SNid,) + —==N'{d2 



-"2 



Irn^T VV2. 



« + '= + ^V2^7M+, + ,, + , + ,V2' 



V2r^^ 1^2/ 



(52) 



where Cbs is given by Eq. (H9ll . d\ and d^ are given by Eq. ( l50ll . N'{x) = dN{x)/dx and 
{}, c;, and k are defined in Eqs. (!38|) - (l40|) . It can be easily proven that the resulting price 
satisfies the so-called put-call parity 



CiS, T, Zo) + Ke-'^ = V{S, T, z^) + S 



where V is the price of the put option which is a derivative contract whose payoff is max(ii' — 
S*, 0) [1]. This relationship guarantees the absence of arbitrage, that is: the practice of taking 
advantage of a price differential between different assets without taking any risk. Finally 
when the initial volatility zq has been averaged out the price has same form as that of 
Eq. (!52|) but using the parameters given by Eq. (145|) . 

IV. SOME RESULTS 

We will now analyze the call price given by Eq. ( I52l) . In Fig. [2] we show the effect of 
changing the parameters of the model on the resulting option price. As we did in Fig. [1] 
each plot slightly modifies only one of the four parameters while the other three are kept 
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FIG. 2: Normalized C/K call price (I52p as a function of the moneyness S/K for T = 20 days 
assuming Aq = 10^'^ and Ai = 10^'^ (cf. Eq. ()19p ) with terms provided by Eqs. (I37p - (l40p when 



zq = 0. We depart from the parameters m = 10 ^ day ' , a = 8 x 10 ^ day , p = —0.4 and 
k = 0.11 day~ ' and slightly modify them in each of these plots. 

constant although taking realistic values. Looking at Fig. [2] we may say that the call price 
is highly sensitive to the normal level m. This is, however, not surprising since an extreme 
sensitivity to volatility -specially around moneyness S/K = 1- is a characteristic feature of 
the classic BS price [l]. As an overall statement we may say that having large values of any 
parameter implies (except for p) a more expensive option. 

We also compare our price, Eq. (!52l) . with the BS price, Cbs^ given by Eq. (H9l) in which 
the volatility is constant. The insets in Fig. [2] represent the difference C — Cbs and we there 
observe a distinct behavior depending on whether we have moneyness smaller than one (in- 
the-money option) or larger than one (out-the-money option). In general, our in-the-money 
calls are cheaper than the BS ones while the out-the-money calls become more expensive 
than the BS ones. This is however true as long as p is negative since otherwise we would 
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FIG. 3: Implied volatility (in yearly units) as a function of the moneyness SjK for T = 20 days 
assuming Aq = 10^'^ and Ai = 10^'^ (cf. Eq. ()19p ) with terms provided by Eqs. (I37p - (l40p when 



zq = 0. We depart from the parameters m = 10 ^ day ' , a = 8 x 10 ^ day , p = —0.4 and 
k = 0.11 day~ ' and slightly modify them in each of these plots. 

have the opposite effect. Let us recall that p should be negative because of the negative 
skewness in the return distribution and also due to the negative return-volatility asymmetric 
correlation (leverage effect), both properties observed in actual markets. The profile of the 
call difference C — Cbs that we have obtained is indeed consistent with the observed one such 
as in Ref. 36| where, although based on a different option pricing method and on a different 
market model as well, the correlation coefficient, p, is studied with special attention. We 
also note that the impact on the call price of changing the reverting force a is much smaller 
than that of changing the vol-of-vol k or the normal level m. 

We now go one step further and study the implied volatility cri. This is the volatility 
that the classic BS formula should adopt if we require that CBsi^'i) = C- We evaluate a^ 
numerically in terms of the moneyness and for an identical set of parameters than those 



15 



to 

s 

T3 





1 1 1 1 1 1 


1 1 11 1 . 


-,U»*^ 


0.8 
0.6 


k=0.11 - 

k=0.05 -- 

-k=0.16 - 


/^ 


■ 


0.4 


y 


/ 


- 


0.2 




1 1 11 1 . 


- 





1 ■ 1 1 ■ ■ 1 


1 . 1 1 1 __l.-i-j^;^ 


0.8 


Ti=0.007 - 
m=0.014 ■■- 
- m=0.01 - 


■0^^^^^' 


0.6 


- 


/ -■ 


0.4 


J 


- 


0.2 


l^ 


1 , . . 1 , . , ' 



to 

CO 
T3 







0.95 1 


1.05 




P= 


-0.4 - 


' ^^^^iii:;^^^^- 


0.8 


:?: 


= 0.0 - y<? 

=0.4 ■■■ /f- 


- 


0.6 


- 


g 


- 


0.4 


_ 


/f' 


_ 






£■' 


■ 


0.2 


- 


/' 


- 




■■•_"--' 


■^ 1 11 1 1 1 . 


■ 





0.95 


1 


1.05 


0.8 


a=0.008 - 

a=0.003 - 

: a=0.04 - 


1 


Z"^"^-- 


0.6 


- 


/ 


- 


0.4 


y 


/ 


- 


0.2 


y 


1 


III.' 



0.95 1 1.05 

moneyness S/K 



0.95 1 1.05 

moneyness S/K 



FIG. 4: Delta hedging ()53p divided by strike X as a function of the moneyness SjK for T = 20 
days assuming Aq = 10^'^ and Ai = 10^^ (cf. Eq. (fT9l) ) with terms provided by Eqs. (f37|) - (H0l ) 



when zq = 0. We depart from the parameters ttt, = 10 ^ day 



^1/2 



a 



8 X lO^^day-^ p = -0.4 



and k = 0.11 day ' and slightly modify them in each of these plots. 



presented in previous figures. We observe that the vol-of-vol, A;, and the long range memory 
parameter, a, have both a rather similar effect although the steepest profile corresponds in 
one case to smaller values (the case of a) while in the other case it corresponds to larger 
values (the case of /c). In contrast, the normal level of volatility m simply shifts the implied 
volatility profile to lower or higher volatility levels it keeps the same form of the profile. The 
smile effect is smirked in one side or in the other depending on the sign of p. The rest of 
the plots studying parameters fc, m^ and a are conditioned to this sign and for all of them 
we take p = —0.4. 

We can also provide an analytical expression for the delta hedging [l|, |27|. This is a 
crucial magnitude since it specifies the number of shares per call to hold in order to remove 
risk of underlying asset price fluctuations (but not volatility fluctuations) from the portfolio. 
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FIG. 5: Call price (I48p and implied volatility (in yearly units) as a function of the moneyness 
SjK for T = 20 days. Left column studies the effects of a non-zero initial volatility assuming 
Aq = 10^^ and Ai = 10~^. Right column shows those caused by changing the constant involved 
in the risk aversion function ()19p when zq = 0. The rest of parameters are m = 10^^ day" ' , 
a = 8 X 10"3 day~\ p = -0.4 and k = 0.11 d&y~^/'^. 



From Eq. (!52|) . it is straightforward to obtain 

dC Ke~^^ 



dS 



S^/rnFT 






2m?T 



H, 



^2 ^ 






(53) 



Figure H] thus provides the same set of plots as those of the previous cases. The long range 
memory parameter a has little effect in the delta hedging although 6 has a higher sensitivity 
to the rest of parameters. Smaller values of the normal level of volatility m make steeper 
the delta hedging profile. The correlation p and the vol-of-vol k have a non-trivial effect 
depending on the moneyness. 
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FIG. 6: Call price as a function of the moneyness S/K for T = 10 days. Points represent the 
empirical call option prices on the Dow Jones Index (DJX) at a precise date (May 2, 2008) and 
with maturity on May 16, 2008. Dashed line takes a call price (j52p fit having fixed the model 
parameters estimated from historical data and with the initial volatility assumed to be the CBOE 
DJIA Volatility Index (VDX) and the current interest rate ratio r = 2%. The curve thus provides 
a fit with proper risk aversion parameters Aq and Ai. 

In Fig. O we analyze the effects on the call price of the initial volatility and the risk 
aversion. Until now we have taken zo = but now we can consider other possible values. As 
expected, we observe in Fig. \5\ that the call becomes more expensive if one takes an initial 
volatility greater than the normal level and cheaper in the opposite case (see also discussion 
in Section [lIIAp . We also look at the risk aversion terms provided by Eq. (TT9|) . Negative 
terms would thus correspond to a more expensive call while positive terms imply having a 
cheaper option since the agent is less risk averse. 

Before concluding this section, we question ourselves whether the risk averse parameters 
Aq and Ai can be in some way or another inferred from empirical data. These parameters 
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give reason of the risk averse perception of the investors and depend on the situation of the 
market at a particular time. Just as an illustrative example, we can look at the European 
option contracts (DJX) on the Dow Jones Industrial Average index traded in the Chicago 
Board Exchange (CBOE) at a given date (2nd of May 2008) and for a particular maturity 
(16th of May 2008, 10 trading days ahead). We next assume that the call price is given by 
Eq. (l52l) with model parameters to be those estimated historical data {m = 10~^ day" ' , 
a = 8 X 10"^ day" , p = —0.4 and k = 0.11 day" ' ), take the US current risk-free interest 
ratio (annual rate r = 2%), and finally consider initial volatility to be the CBOE Volatility 
Index (VDX) on the 2nd of May 2008 (annual rate o"o = 16.55% from which we get yo = 
ln(o"o/m)). The estimation on o"o can be more or less sophisticated |23| but it is rather usual 
to assume the VDX since this is designed to reflect investors' consensus on current volatility 
level. Figure [6] therefore presents a rather satisfactory flt on the empirical option contracts 
over different moneyness S/K by solely modifying the risk parameters Aq and Ai. The error 
in the fltting is rather small being of the same order or even smaller than the tick size of 
the traded option. In this way and thanks to the fact that the parameters of the model are 
obtained from historical Dow Jones data, we can obtain the implied risk aversion parameters 
which typically fluctuates in time as investors' perception changes. 

V. CONCLUDING REMARKS 

The main goal of this paper has been to study the effects on option pricing of several 
well-known properties of flnancial markets. These properties include the long-range memory 
of the volatility, the short-range memory of the leverage effect, the negative skewness and 
the kurtosis. The analysis is based on a market model that satisfles these properties which, 
in turn, can be easily identifled through the parameters of the model. In this way we provide 
a different and more complete analysis on option pricing than that we had presented some 



time ago in which the effects of non-idea 
relaxation were taken into account 37 



market conditions such as fat tails and a small 

Si. 



We have derived an approximated European call option prices when the volatility of 
the underlying price is random and it is described by the exponential Ornstein-Uhlenbeck 
process. The solution has been obtained by an approximation procedure based on a partial 
expansion of the characteristic function under the risk-neutral pricing measure. 
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The call price obtained is valid for a range of parameters different than those of a previous 
study on the subject [8]. In that work Fouque et el assumed that the reversion toward the 
normal level of volatility is fast. In other words, the parameter a is large and the character- 
istic time scale for reversion, l/a, is of the order of few days. Fouque et el also considered 
that (3"^ = k'^/2a ~ 1 which implies that the return- volatility asymmetric correlation {i.e., 
the leverage effect) should have a characteristic time comparable to that of the volatility au- 



tocorrelation. Mostly based on the Dow Jones daily index data [13, l23j, we have considered 
a rather different situation where the fast parameter is not the reverting force a but the 
vol-of-vol k. In this way, we have singled out these two market memories thus allowing for 
a leverage effect during a time-lapse of few weeks and a persistent volatility autocorrelation 
larger than one year. These properties are consistent with empirical observations on the 
Dow- Jones index 12 1. 



Under these circumstances we have constructed an approximate option price where risk 
aversion is assumed to be a linear function of the logarithm of the volatility. This approx- 
imation contains corrections in the variance, the skewness and the kurtosis, all of these 
corrections in terms of Hermite polynomials. This constitutes a tangible step forward with 
respect to other approaches which use the Heston stochastic volatility model [l5|, [l9(] but 
only consider zeroth-order corrections. Our a ppr oach to the martingale measure albeit be- 



ing more complete than those taken in Refs. |15l . Il9l | it is still able to provide an analytical 
expression for the call price and the subsequent Greeks. 

Summarizing, we have studied the call price and its implied volatility and observed that 
the correlation p between the two Wiener input noises plays a crucial role. The behavior of 
the call can greatly change depending on the sign of p which confirms the findings of the 
previous work of Pochart and Bouchaud [3&\ and many others. We have therefore focused 
on a negative value of p that is consistent with empirical observations of the leverage effect. 
Keeping p constant, moderate values of risk aversion, and a maturity time of the order of 
few weeks, we have also observed that the vol-of-vol, k, and the normal level of volatility, 
m, have both an important impact to the option price. This appears in clear contrast with 
the very little effect of not having a reliable estimation of the rate of mean-reversion of the 
volatility quantified by the parameter a. 
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APPENDIX A: APPROXIMATE SOLUTION OF THE CHARACTERISTIC 
FUNCTION 

We start from Eqs. fl33l) - fl34l) and look for an approximate expression of the joint distri- 
bution (p{uji, UJ2, t') valid for large values of A. We also note that the marginal characteristic 
function of the (martingale) return can be obtained from the joint characteristic function 
by setting 002 = 0. Therefore, we will look for a solution to the problem (!33|) -( l34|) that for 
small values of uj2 takes the form: 

^{ui, UJ2, t') = exp {- [A{ui, t')ujl + B{ui, t')^2 + C{uji,t') + 0{ul l/A^)] } . (Al) 

Substituting this into Eq. (133|1 yields 

AujI + BUJ2 + C = -2vujIA - VUJ2B + —ujI - tt^ + ^ + -ujI 1 + 4— ^v4 + —B 



2 ' 2X m^X 2 ' \ X X 



+XpuJiUJ2[l + ^A+jB)+0{l/X'), (A2) 



where the dot denotes a time derivative. 

Collecting the quadratic terms in ti;2, we get 



\2 

A = -2{u~ipuJi)A + —. 



The solution to this equation with the initial condition A{t' = 0) = is 

A^ 

47 (wi) 



A{uj,,t') = ^-- [1 - e-2T(^i)*'l (A3) 



where 

7(c<Ji) = z/ — ipui. (A4) 

We now take the linear terms in 002'- 

2icj2 

B = -75 H -^A + Xpuji. 

X 

The solution to this equation with the initial condition B{t' = 0) = —Wq is 

B[u,,t') = -..oe-*' + ^ (1 - e-*') + ^ (l " e-*')' (^5) 

where 7 = 7(ci;i) given by Eq. flA4p . Finally, the terms independent of UJ2 yield 
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with C(t' = 0) = 0. Therefore 
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27 V / 7 

We have thus obtained all the terms in the two-dimensional characteristic function deter- 
mined by Eq. (155]) . 

However, that we only need to know the marginal characteristic function of the return 
X{t) which implies that we only have to solve the equation when uj2 = and assume 
Ui = uj/\. From Eq. (jAip we have the approximation 



<^{uj/X, 0, t') = exp {-C{uj/X, t')} 



(A7) 



where 

C{uj/X,t') 
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2A472 
Note that there is an extra parameter involved : 

'j{uji = uj/X) = V 
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■7*' 
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that depends on A and it leads us to write CiiojX^ t') in a somewhat more compact form (cf. 
Eqs. f lA4p ). Indeed, the first term we have to reconsider is 
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the second one reads 



A37 



t'--(l- e-^*' 
7 ^ 

1 



A3i/2 
1 

A3^ 



ut' 



1 - e--') + ^ 
^ Az/ 



(z/t'(l + 



e--*') - 2(1 - e-"*' 



+ ^) .t' - f 1 + ^) a(t') + ^6(0 
Az/ / V Az/ / Az/ 



+ 0(1/A^) 
+ 0(1/A^); (AlO) 



22 



while the third one is 
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(All) 



where 

a(t') = l-e-"*', 6(t') = ^^t'e""*'. 

All these expressions serve us to study the terms included in C{uj^ t') given by Eq. (1A8I) up 
to order 1/A^. We sum up the contributions flA9p - flAlip and obtain 
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We finally rearrange this expression taking into account the order of u. The final result is 
shown in Eq. ( l36l) of the main text. 



APPENDIX B: DERIVATION OF THE EUROPEAN CALL OPTION 



We perform the average given by Eq. psj) . Due to the fact that we have four contributions 
in Eq. (H2!) . we will also obtain four terms for the option price. We decompose them as follows 



C{S,T,zo)=Cbs{S,T)+^{T,zo)Co{S,T)+pc;{T,zo)Ci{S,T) + 



^{T) + -{}{T,z,f 



C2(5,T), 

(Bl) 
where first term Cbs corresponds to the Black-Scholes price {i.e., when underlying process 
has a constant volatility given by m). The following terms -containing several corrections 
in the volatility, the skewness and kurtosis- can be easily derived if one considers 



e-"^ Hn{-a) 



da^- 



(B2) 



where Hn{-) are the Hermite poljTiomials. The first term due to a non constant volatility is 
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and taking into account Eq. (lB2p it becomes 
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which, after some manipulations, finally reads 

Co{S,T) = SN{d,) + ^==N\d2), 

where N'{x) = dN{x)/dx and di and d2 are defined in Eq. fl50l) . 
The second term of our calculation is 
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X — fl 



Again taking into account Eq. (1B2I) . we have 



exp 
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i/i (t/2/V2) - 1 



(B4) 



which, after some simple algebra, yields 

C,{S,T) = SN{d,)-^==N'{d,) -j== 

Note that correlation p between the two Brownian noise sources determines the sign and the 
strength of this term. Obviously if there is no correlation this term disappears (cf. Eq. (JBll) ). 
The third and last piece of our option price reads 
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and, after using Eq. (]B2p . we get 
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which yields 



C,{S,T) = SN{d,) + ^==N\d2] 

ym'^l 



We sum up all contributions given by Eqs. fIBSp - flBSP and plugging them into Eq. fIBip we 
finally obtain Eq. (1521) of the main text. 
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